Particle corrélations in a fermi superfluid 
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We discuss corrélations between particles of différent momentum in a superfluid fermi gas, ac- 
cessible through noise measurements of absorption images of the expanded gas. We include two 
éléments missing from the simplest treatment, based on the BCS wavefunction: the explicit use of 
a conserving approximation satisfying particle number conservation, and the inclusion of the con- 
tribution from Cooper pairs at finite momentum. We expect the latter to be a significant issue in 
the strongly correlated state emerging in the BCS-BEC crossover. 
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The récent observation of résonance pairing in gases 
of fermionic atoms in the BCS-BEC crossover régime [ï| 
has been followed by an explosion in expérimental and 
theoretical activity aimed at understanding this canoni- 
cal problem of many-body physics. The measurement of 
the paring gap 0] and most recently the observation of 
superfluid flow [3[ have dramatically confirmed the prés- 
ence of superfluid condensâtes in thèse Systems. 

One experiment that is yet lacking is a direct mea- 
surement of order, that is, a démonstration of conden- 
sation of pairs into the a zéro momentum state. In the 
case of Bose-Einstein condensation of bosonic atoms (or 
molécules), the distribution n(k) of the number of parti- 
cles in each momentum state is readily available by time- 
of-flight imaging. The signature of BEC is a sharp peak 
in this distribution at zéro momentum. By contrast, n(k) 
in a fermionic superfluid is a smooth function that re- 
sembles the corresponding distribution in a normal fermi 
gas at finite température - see the récent experiment [ij . 
Thus such a one-particle measure is of little use. 

A ingenious resolution to this problem was the sugges- 
tion in Réf. [j| to measure noise corrélations in the time- 
of-Hight images. The density-density corrélation function 
of the gas following expansion can be interpreted as the 
(column-integrated) correlator of n(k) at two différent 
momenta in a gas prior to release from the trap. Then 
the condensation of pairs at zéro momentum corresponds 
to a peak in the correlator 



to give 
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at opposite momenta ki = — k 2 . We dénote the two 
atomic species between which pairing occurs as s =î, |, 
and assume species-selective imaging. In Réf. Q this 
corrélation function was evaluated (at zéro température) 
using the usual BCS ground state ansatz 
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where E p = y^ 2 + |A| 2 , and £ p = p 2 /2m - p. The 

first experiments demonstrating the possibility of such 
measurements have now been performed @, 0] • Since a 
measurement of Cfj(ki,k2) in the BCS-BEC crossover 
may be attempted soon, and given the utility of such a 
measurement in constraining many-body théories of the 
strongly interacting gas, our purpose in this paper is to 
develop the theory of noise corrélations in the following 
two directions: 

1. Conserving approximation. We note that Eq. ([3]), 
together with the correlator C||(ki,k2), which has the 
same form but with a delta- function #ki,k2j ca n be used 
to calculate the density structure factor 5*(q) = (pqp-q) 
of the System at zéro momentum. The resuit (|ÏÏJ gives 
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(in the BCS limit, where v is the density of states per 
atomic species at the fermi energy) However, it is well 
known [| (see for a récent discussion in the contcxt of 
fermi gases) that the structure factor of a quantum fluid 
with finite compressibility vanishes at zéro momentum at 
zéro température. The incorrect resuit (j4|) first appeared 
in [I0( , where it was used to conclude that the low energy 
modes of a neutral superconductor would have dispersion 
uj ~ q 2 , rather than the linearly dispersing Bogoliubov- 
Anderson mode. Since (U) arises entirely from the corré- 
lations found in Réf. one may worry that the resuit 
([3]) will be changed in some essential way. 

Note that this problem is unrelated to the non- 
conserving form of Eq. ((2]). It is easy to see that the 
projected form 
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lacks the long-ranged corrélations in position space re- 
quired to make the structure factor vanish }16j . 

The resolution to this problem, as described in the 
book is to use a conserving approximation. In the 
language of diagrams, we must include vertex corrections 
rcquired by particle number conservation. Our first goal 
in this paper is to explain the effect thèse have on the 
correlator Cj j(ki,k 2 ). The good news is that the delta- 
function structure of (0]) remains intact, but the corréla- 
tions function acquires a smooth background 



Cu(ki,k 2 



C î B | cs (k 1 ,k 2 )+C(k 1 ,k 2 ), (5) 



where C(ki,k 2 ) dépends only on the magnitude of the 
two momenta, and is smooth on the scale of the gap. 
This term cancels the first one when we compute the 
structure factor as in Eq. 2] It is natural to interpret 
C(ki,k 2 ) as arising from corrélations between différent 
pairs so that, as in a single-species BEC, the vanishing of 
the structure factor at zéro wavenumber is due to pair- 
wise corrélations between bosons (i.e. the Bogoliubov 
approximation). Such an interprétation implies that a 
measurement of /our-particle corrélations would reveal a 
peak at ki + k 2 + k3 + k4 = 

2. Cooper pairs at non-zero momentum. Another well- 
known conséquence of including corrélations between 
bosons is that the distribution function n(q) of the num- 
ber of bosons with given momentum acquires non-trivial 
structure (leading to condensate depletion). Thus the 
true corrélation function C||(ki,k 2 ) will likcly have a 
strong angular dependence outside of the main peak at 
ki = — k 2 duc to the présence of such pairs, particularly 
in the crossover région where we expect déviations from 
the BCS picture to be most extrême. In fact, part of this 
dependence turns out to be universal and has the form 
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where N p = k F /(6ir 2 ) x Vol, with Uf the fermi wavevec- 
tor, is the number of pairs i.e the number of atoms of each 
species, c s is the sound velocity, and £ is the Ginzburg- 
Landau cohérence length that interpolâtes between the 
familiar cohérence length /cp/mA of the BCS superfluid 
and the healing length l/4mc s of the BEC [12]. The 
meaning of this contribution is very simple (see Fig. ^jp). 
Outside of the brackets we have the square of the usual 
Cooper pair wavefunction 0(k) = A/Ek, describing the 
internai structure of the pairs. Inside the bracket we 
have a new contribution, which should be interpreted 
as the the momentum distribution of uncondensed pairs 
n p (q) — (2m)c s /2q. This form is precisely that expected 
m simple BEC [3, El, and can be deduced very gener- 
ally from the uncertainty principle [Ï5j . 

The contribution of pairs of finite momentum is al- 
ready one of the issues in the experiment [|| , which mea- 
sured the corrélations of dissociated molécules. The re- 
suit ^ shows how superfluid order manifests itself in 




FIG. 1: Diagrams corresponding to a) the naive approxima- 
tion Eq. (|3]) and b) the vertex correction required for a con- 
serving approximation. 



this contribution. Clearly this term is small in the ex- 
trême BEC limit as c s — ► 0. In the BCS limit c s goes to 
the finite value vf/V3, where Vf is the Fermi velocity, 
but the région of validity of ^ becomes very small as 
£ — > oo. Sincc £ attains a minimum value of order the 
fermi wavelength in the crossover this is a promising place 
to look for the contribution of finite momentum pairs. It 
would be very interesting to experimentally confirm thèse 
strong déviations from the prédiction of the simple BCS 
wavefunction in the crossover state. 

We now turn to the dérivation of thèse results. Let 
us first repeat the calculation of Réf. la in the Green's 
function language to establish our formalism. In terms of 
fermion création and annihilation operators the number 
of particles of each species with momentum k is n s (k) = 
if>l(\i)ip s (\s.). To calculate the corrélation function fi} we 
introduce the Gor'kov Green's functions (real time, zéro 
température) 

G(ki, h; k 2 , t 2 ) = -*(Ttf (ki, h) <g> *t(ka,f 2 )) 



in terms of the Nambu spinor 
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Evaluating the corrélation function |T]) using Wick's 
theorem gives 

C n (ki,k 2 ) = -Gi 2 (ki,0;-k 2 ,0)G 2 i(-k 2 ,0;ki,0) 

The diagram corresponding to this contribution is 
just the bubble (see Fig. [T]). The explicit form of 
the Green's function is G(ki,£i;k 2) £ 2 ) = {2-k)6{e\ — 
£2)4^2 Go(ki,£i), with 
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Straightforward calculation yields the resuit ^ of Réf. 5. 
The delta-function peak at k x = — k 2 is the signature of 
condensation of Cooper pairs with opposing momenta. 
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We could also study the corrélations of the same species: 



C TT = Gu(ki, +; k 2 , 0)Gn(k 2 , 0; k x , +) = 5 klM 



4^, 



(we have to pay attention to the point splitting in the 
time direction to get the right answer). The origin of 
this term is simply the binomial variance n(k) (1 — n(k)) 
arising from having < n(k) < 1 due to the pairing, so 
this part resembles the corrélations at finite température 
in a normal fermi gas. 

As explained earlier, the resuit ([3]) cannot be ail there 
is, as it gives the wrong value of 5(0). In gênerai, 
the approximation for a two particle quantity such as 
C||(ki,k2) must be consistent with that used for the 
self-energy. The self-energy approximation in question 
is just the usual BCS self-consistency équation. We 
now describe the resulting approximation for C|j.(ki, ka). 
Digrammatically, vertex corrections (see Fig. [T]) are re- 
quired to ensure particle number conservation, i.e. that 
the vertex function satishes the generalized Ward iden- 

tii.v [m. 

Consider the case of the model interaction 
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We can represent this operator in Nambu form in the 
following way 



- (* (g) *t| T3 (g, T3 |$t (g, ijr) , (8) 



(some dclta-function terms arising from anticommutation 
drop out). The final form is for later convenienec. The 
simplest approximation for the self-energy corresponding 
to this interaction is 



E(p,e) = iA 
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Note that this includes the Hartree and Fock channels 
as well as the Cooper channel. The Hartree part is the 
second term of ©; the Fock part is included in the first 
part. The Hartree- Fock potential appears as the coeffi- 
cient of T3 in the Cooper part is off-diagonal, and 
reproduces the usual self-consistent équation. In the fol- 
lowing calculation, we will work in the weakly coupled 
(BCS) limit. For now we drop the Hartree part (second 
term), and will explicitly demonstrate shortly why the 
Cooper part is the only one that needs to be accounted 
for in the weakly coupled problem. 



This first term of Eq. ([9]) implies that the vertex func- 
tion should be taken to be the géométrie séries of bubbles 
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The matrix structure of (|10p is inherited from ((H}. A(q, to) 
is a 4 x 4 matrix. We dénote the restriction of Aa^^i to 
the subspace a ^ b', a' ^ b as A c , corresponding to the 
Cooper channel. The restriction <3E> c (g,cj) is afflicted with 
the usual log-di vergence. To handle this, we use the gap 
équation 
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(y is the fermi surface density of states per spin degree 
of freedom) to write 



v \ de 



i 



VA2+F (4 (A 2 + i 2 ) - lo 2 ) 

-A 2 + lo 2 /2 -A 2 
-A 2 -A 2 +lu 2 /2. 



Outside of the Cooper subspace, 5>(0,u/) is completely 
convergent, so that when we taken the limit A — * as 
the cut-off goes to infinity (limit of a delta-function po- 

tential), we can treat 1/A — <ï> 
subspace. 

The final form of the vertex function is 



as zéro outside this 
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Note that the vertex matrix has one eigenvalue A^O, lu) — 
A2(0,u;) that goes like 1/lo 2 as ui — > 0. This is due to the 
phase mode. The other corresponds to the amplitude 
mode. 

It is now straightforward to find the vertex correction 
to Cfj.(ki,ka). We have to compute 
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Y, [ dtidhfâfa, 0)* x (ki, 0) (*t(k, fi) ® *(k, ii)|A(0, ii - t 2 ) |*(k', fa) ® *t(k,' ta)) tf 2 (k 2 , 0)* 2 (k 2 , 0)J 
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The vertex correction shown in Fig. Q] corresponds to pairing the V&i's with one side of the vertex operator, and the 
with the other side. The other diagram, that corresponds to the 'crossed' pairings, will be discussed later. Taking 
into account that only the 2 — 3 subspace of the vertex operator is non-zero gives 
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G u (£1 + w/2, ki) G n {ex - w/2, ki) G j2 (e 2 + w/2, k 2 ) G 2j (e 2 - w/2, k 2 ) A CÎJ (0, w). (12) 



(1 = 2, 2 = 1). Performing the £i i2 intégrais yields the contribution of the vertex correction to C^(ki,k 2 ) (and also 
toC n (ki,k 2 )) 
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FIG. 2: a) Diagram describing contribution from finite mo- 
mentum Cooper pairs, b) Schematic of the disintegration 
process 



cancels (|4"|). Turning now to the contribution of finite 
momentum Cooper pairs, we will see that the diagram 
in Fig. [2k is the leading contribution to C|j_(ki — > — k 2 ). 
This diagram simply corresponds to the 'crossed' pairing 
in (fTTj) . We can ignore vertex corrections, as we know 
that they give rise only to a background independent 
of the angle between ki and k 2 , so in particular don't 
contribute to any singular contribution. The resulting 
expression is in gênerai rather lengthy. Note that we 
now require the vertex function at finite q = ki + k 2 . At 
low q, however, the final f2 intégral is dominated by the 
singular residue of the phase mode Ai(q, w) — A 2 (q, u>) 



The resuit does not dépend on the relative angle of ki 
and k 2 . Now we can check the basic idea: that this con- 
tribution cancels the naive one calculated in Eq.[3]to give 
a vanishing structure factor at zéro wavenumber. First 
note that after £ intégration only the Ai — A 2 (phase 
mode) part of the vertex function contributes. Perform- 
ing the £ intégrais gives. 
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With the required factor of 4 for the spin sum (TTJJ exactly 



Res [Ai(q,w) - A 2 (q,w] w=± = ± 



2mc s A 2 
n p q 



where c s = Vf/V3 and n p — kp/6ir 2 is the density of 
pairs. This singular behavior gives rise to the second 
term in ©. Though we worked in the BCS limit, we 
believe this form is valid throughout the crossover, with 
the appropriate sound velocity c s . Note that because we 
worked at weak coupling, no depletion of the conden- 
sate is accounted for. With depletion, both terms in (J6]) 
should be scaled by the condensate fraction. 

The author would like to thank Victor Gurarie, Léo 
Radzihovsky and Dan Sheehy for discussions, and for 
their hospitality during a visit to University of Colorado, 
where this work was completed. 



[1] C. A. Régal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. [2] C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, 
92, 040403 (2004). S. Jochim, J. H. Denschlag, and R. Grimm, Science 305, 



■5 



1128 (2004). 

[3] M. W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, C. H. 

Schunck, and W. Ketterle, Nature 45, 1047 (2005). 
[4] C. A. Régal, M. Greiner, S. Giorgini, M. Holland, and 

D. S. Jin, cond-mat/0507316 (2005). 
[5] E. Altman, E. Dernier, and M. D. Lukin, Phys. Rev. A 

70, 013603 (2004). 
[6] M. Greiner, C. A. Régal, J. T. Stewart, and D. S. Jin, 

Phys. Rev. Lett. 94, 110401 (2005). 
[7] S. Fôlling, F. Gerbier, A. Widera, O. Mandel, T. Gericke, 

and I. Bloch, Nature 45, 1047 (2005). 
[8] P. J. Price, Phys. Rev. 94, 257 (1954). 
[9] H. P. Buchler, P. Zoller, and W. Zwerger, Phys. Rev. 

Lett. 93, 080401 (2004). 
[10] P. W. Anderson, Phys. Rev. 110, 827 (1958). 



[11] J. R. Schreifïer, Theory of Superconductivity (Perseus 

Books, Reading, MA, 1964). 
[12] J. R. Engelbrecht, M. Randeria, and C. A. R. Sa de Melo, 

Phys. Rev. B 55, 15153 (1997). 
[13] J. Gavoret and P. Nozieres, Ann. Phys. (NY) 28, 349 

(1964). 

[14] G. V. Chester and L. Reatto, Phys. Lett 22, 276 (1966). 
[15] L. Pitaevskii and S. Stringari, J. Low Temp. Phys. 85, 
377 (1991). 

[16] We should point out that the finiteness of S'(q) at q — > 
is not inconsistent with the value exactly at q = being 
zéro, as it will be for a number-conserving wavefunction. 
The requirement S(q — > 0) — ► is a stronger constraint, 
not met in gênerai by such wavefunctions. 



